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We study photon production from the hidden local symme- 
try approach that includes n, p, and ai mesons and compute 
the corresponding photon emission rates from a hadronic gas 
in thermal equilibrium. Together with experimental radiative 
decay widths of the background, these rates are used in a rel- 
ativistic transport model to calculate single photon spectra 
in heavy-ion collisions at SPS energies. We then employ this 
effective theory to test three scenarios for the chiral phase 
transition in high-temperature nuclear matter including de- 
creasing vector meson masses. Although all calculations re- 
spect the upper bound set by the WA80 Collaboration, we 
find the scenarios could be distinguished with more detailed 
data. 
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I. INTRODUCTION 

Experiments due to be performed over the next few 
decades are aimed at achieving a quark-gluon plasma, 
thereby hopefully granting the ultimate confirmation of 
QCD as the theory of strong interactions and also settling 
the debate over the phenomenology of hadronic matter 
at high density. However, data currently available from 
experiments carried out at the CERN SPS in the form 
of dilepton and photon spectra offer information relevant 
to many outstanding issues, such as the enhancement 
of low-mass dileptons in central heavy-ion collisions re- 
ported by the CERES and HELIOS collaborations §,|, 
anomalous J/ip suppression in central Pb-|-Pb collisions 
by the NA50 collaboration Q, and finally the limit on 
single photon spectra from central heavy ion collisions 
set by the WA80 collaboration ||] . 

There has been intense theoretical activity in the past 
few years to predict the behavior of hadronic matter at 
high density as chiral symmetry restoration and decon- 
finement are approached. However, the general physical 
picture of QCD phase transitions is still not agreed upon. 
Lattice QCD lacks the techniques to surmount the finite- 
density problem and other models differ in opinion to the 
important effects. In particular, there is nothing to pro- 
hibit chiral symmetry from being restored via the Georgi 
vector limit [|j . 

In this limit, the p becomes massless. Its longitudinal 
polarization must disappear, transforming into a scalar- 
isovGctor particle identified as the chiral partner of the 
pion. The ai is the chiral partner of the p and so also 



must become massless. The interesting question is the 
nature of the mass decrease of the vector mesons and 
whether heavy ion collisions will be able to distinguish 
between them. The decrease in general is an effect still 
in contention although it is able to explain the dilepton 
spectrum from CERES within a transport code calcula- 
tion §. 

Based on the restoration of scale invariance of QCD for 
low momentum scales. Brown and Rho 1^ suggested that 
the mass of non-strange vector mesons should decrease in 
dense matter, together with the chiral condensate. There 
are a number of arguments that emerged that are in favor 
of Brown-Rho scaling ||9[]. However, finite-temperature 
calculations based on effective Lagrangians |lQ,Ol and 
some lattice results |1^ find that the p mass does not 
vanish in the chiral linrit. 

The aim of this paper is to compare the predictions 
of decreasing in-medium vector meson masses (whether 
relying on the Georgi vector limit or not) with the 
WA80 photon spectra. The most important processes 
are tttt — > p^ and irp — > 717 as well as the decay p — *■ tttt"/ 
as determined in Ref. |13| ]. The ai resonance is also im- 
portant for photon rates, as discussed in |14-|l9], and 
was treated incompletely in the transport calculations 
of Ref. |l^. Furthermore, different scenarios for the den- 
sity dependence of matter such as the Georgi vector limit 
have not been addressed in the context of photon rates 
in the past. We therefore use the extended hidden local 
symmetry Lagrangian to make easy tests of different sce- 
narios through the adjustment of the parameters. The 
Georgi vector limit is possible in the hidden gauge de- 
scription as described in the next section. The massive 
Yang-Mills Lagrangian, as used in ||l^, is to some ex- 
tent a special case of this Lagrangian, corresponding to 
a specific parameter choice and a different gauge-fixing 
scheme. 

In general, the effective parameters - masses and cou- 
plings - may be seen as a result of two distinct steps. 
First, the low-energy limit of QCD defines the physics of 
hadrons in free space due to chiral symmetry breaking 
and confinement. Second, the physics of hadrons at high 
density and temperature is obtained through many-body 
effects with the above mentioned low-energy interaction. 

One might argue to what extent this dichotonry is a 
natural thing. Whether it is suggested by the very ex- 
istence of hadrons, or is a result of the historic devel- 
opment of strong interaction physics. There is no good 
reason to exclude medium effects from either of these 
steps. For instance, finite-temperature lattice QCD cal- 



culations, in principle complete, consider ensembles of a 
very small number of particles. Temperature effects at 
this level are mainly due to modification of the QCD vac- 
uum, and may be compared to internal structure modifi- 
cations of hadrons, rather than many-body effects. Tra- 
ditional many-body calculations, on the other hand, use 
the same parameters in their basic Lagrangian at all tem- 
peratures and densities. Each of these approaches taken 
separately might miss part of the picture. 

In this paper, we do not attempt to predict medium 
modifications of the vector meson masses. Instead, we 
assume all medium effects are induced by Brown-Rho 
scaling of the vector meson masses and the pion decay 
constant through three different scenarios as presented 
below. We therefore can work both in free space and 
at high temperature and density by adjusting the cou- 
plings within the hidden gauge theory framework as will 
be shown in the next section. 

The remainder of the paper is structured as follows. 
In Section II we briefly review the hidden local symme- 
try (HLS) formalism extended to both vector and axial- 
vector mesons and its fundamental results. Then we 
propose a scheme of parameter assignment which leaves 
only the vector meson masses as external parameters. In 
Section III we present our main analytical results, the 
p — > 7r7r7 decay width and cross sections for tttt — > p7 
as well as irp —^ n^, and show the thermal equilibrium 
photon rates derived from them. In Section IV we dis- 
cuss different ways to implement Brown-Rho scaling and 
their effect on the results of the previous section. In Sec- 
tion V we describe the transport model we use and then 
present the predictions for photon rates from S-|-Au col- 
lisions with the different implementations of Brown-Rho 
scaling. Finally, Section VI is devoted to a summary and 
conclusions. 



II. HIDDEN LOCAL SYMMETRY MODEL 

We start with a brief review of the extended HLS 
model. The reader is referred to |18| for further details. 
We would like to include the dynamics of both the p 
and ai mesons and at the same time have a mechanism 
which allows for the Georgi vector limit H. This has 
been done for the p meson alone using an SU{2) hid- 
den local symmetry 113]. In that case the Lagrangian 
depends on a single parameter a with the KSRF rela- 
tion 
the t 



1^ gp — 2/^(7p7r7r being a general consequence of 
tieory. Imposing a — 2 gives the universality of 
the p-couplings, a second KSRF relation m^ = 2g'^^^f^, 
and p dominance which are all phenomenologically mo- 
tivated. The Georgi vector limit is approached as a ^ 1. 
In this limit, the ordering of the hidden local symme- 
try Lagrangian in terms of the number of derivatives can 
be vindicated since both the tt and p are light. Away 
from the vector limit this ordering of terms can still be 
considered valid. 



In order to include the ai, we need to extend the lo- 
cal symmetry to SU{2)l x 5'f7(2)ij. This increases the 
number of terms in the most general Lagrangian to give 



C = a Cy + h La + c Cm + d C-n + C 



kin. 



(1) 



with the kinetic terms for the vector V^ and axial-vector 
A^ fields included in the last term. In addition we can 
add a global symmetry SU{2) x SU{2) to describe ex- 
ternal vector Vp (including the photon) and axial-vector 
An fields. 

Defining the conventional non-linear sigma field for 
the pion in terms of three SU{2) fields: U{x) — 
i\{x)^M{x)£,B{x), we have ||l|l 
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D,.Cr,l = d^in.L - ig{V ± A)niR,L + ieS,R,L{V ± A)^ 

The choice d= \ — h + 6^/(6 -|- c) ensures the pion kinetic 
term in C has unit coefficient. Fixing the gauge ^m = 1 
and £,L — ^R — 6xp(i7r//7r), there is ir-ai mixing present 
in Ca which can be eliminated by the shift 

A^ -> Af, + —— -j^df^TT. 
b + c Ug 

The gauge choice for A and V could introduce two 
other unphysical scalar fields (ct and p) that parameter- 
ize the decomposition of U{x) above. This also leads to 
mass generation for these new fields through the usual 
Stueckleberg construction. However, choosing the uni- 
tary gauge rids the Lagrangian of the unphysical fields 
at merely the expense of making the vector-field propa- 
gators non-transverse. Also, the vector part reduces to 
the CCWZ effective Lagrangian in this gauge [|| . We will 
stick to this gauge throughout for convenience. 

Unfortunately, upon adding the ai, the first KSRF re- 
lation develops a strong momentum dependence and the 
decay width Ta-^^TT'y vanishes. In addition, the width 
Tai^pTT is half of its experimental value. (For a full 
discussion of these points see |l^.) These unattractive 
features can be rectified by including higher derivative 
terms. The added part is chosen to be pSl 
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strength for the field combination (F ± A)^. For a = 2 
the first KSRF relation gp — 2f^gpTr^ and vector domi- 
nance are regained. The oi decays take on the reasonable 
values Ta.^pn = 360 MeV and Ta.^jn = 320 keV. 

Saturating the two Weinberg sum rules with the nar- 
row width approximation [pO[ gives the two relation^j 



9p = 5ai 



m„ 



"Y^P 



(2) 



which are satisfied for b = a and c — a/{a— 1) in eq. M). 
This leaves only the one parameter a to govern the evo- 
lution of the Lagrangian towards the Georgi vector limit. 
For a = 2, the familiar relation m^^ = 2?7i^ is repro- 
duced. This gives rua^ = 1090 MeV whereas the ex- 
perimental value is ttt-qi — 1230 MeV. This discrepancy 
suggests that a = 2 is not the proper value in free space 
and it should be closer to a = 1.64. 

Setting a — 1.64 changes the other physical observables 
as well. First of all, the jmr vertex is not entirely vector 
dominated, resulting in a non-zero g^Tr-n- = 0.18e direct 
coupling to the photon. (The total contribution including 
the p-j mixing is still, of course, equal to e.) The p- 
coupling universality is modified to become ^p^rTr = 0.96(7, 
still within reason. In fact, the two KSRF relations can 
be written for general a as 



gpTTTT = 2"-9(l + ^r) 



a(l + er)' 






f 



2, (r2)- = 0.39 
0.44 ±0.03 fm^. 



with £_r = 2r(l — 2r) and r defined in eq. (||). Using 
/^ = 92.4 MeV and 5p^^/47r ~ 3.0 [|l) improves the 
second KSRF relation prediction for the p mass from 
mp = 800 MeV (a = 2) to 760 MeV (a = 1.64). Notice 
that Georgi's relation m^ = ig^^^f^ is reproduced for 
a = 1. The pion vector radius for a 
fm^, is slightly below the data (r^)y 
For general a the value is 

, _ 3(1 +g.) 
^^^~ 9'fS 

which for a = 1.64 results in 0.46 fm^ in good agreement 
with the data. The pion polarizability can then be de- 
termined from the Das-Mathur-Okubo relation [g2| and 
gives values consistent with data for both values of a. 

The decays of the oi meson also are modified. The 
total width is dominated by Fa^^^p^ which changes from 
360 MeV to 430 MeV with a = 1.64, stiU about equal to 
the experimental value of 400 MeV. The surprise is the 
radiative decay Tai^-yTv- It changes from 320 keV to 50 
keV as given by the formula 



^Adding widths to the resonances changes the quantitative 
result for the parameter a by less than a percent. 



^ ai — ^77T [ 



(3r-l)^|p|=^ 



This is very different from the experimental value 
630 keV also obtained by vector meson dominance. This 
is due to the contribution and interference from the graph 
with a direct photon. This is responsible for pulling the 
width down to 320 keV with a — 2 {r = 1/2) and almost 
vanishing for a ~ 1.64 (r ~ 1/3). 

Therefore, not assuming total vector dominance in 
nature actually improves most quantities towards their 
measured values by about 10% with the exception of the 
decay width Fq^^^^ which nearly vanishes to O(e^). 



III. ANALYTICAL RESULTS 

Although we will also consider the effects of the cu, rj, 
and rj' mesons, these particles are long enough lived to 
mostly decay after freeze-out and therefore we can take 
the experimental values for their partial widths. The 
only decay abundant enough to consider is then p — > mrj. 
For the photon production from two-body collisions, only 
TTTT — > p7 and irp —f nj are important enough to consider 
for the temperature region relevant to heavy-ion colli- 
sions at SPS energies Q. These three processes are 
all related to the same matrix element quoted in Ap- 
pendix A. Electromagnetic gauge invariance can easily 
be checked on our expressions. The contributing graphs 
are shown in Fig. ^ 
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FIG. 1. The general diagrams for the matrix element in- 
volving prnvy. 



First focusing on the decay p 
space for a = 2 (1.64) is 



nn^ 



the result in free 



r(p° -^ TT+n--f) = 1.6 (1.8) MeV 

Both values are in reasonable agreement with the exper- 
imental data 1.5 ± 0.3 MeV [El|. In the same model, 
we can calculate the radiative decay width of charged 
rho mesons, and we obtain r(p -^ ttttj) = 0.88 (0.98) 
MeV after isospin averaging, showing the p^ together 
give slightly less contribution than the p° alone. 

The cross sections for tttt — > p"f and np — > t:^ can be 
calculated in the same way, and the results are shown in 



Fig. H and || for three values of a. A momentum depen- 
dent width for the ai mesons as calculated in the HLS 
Lagrangian was included. It is useful to compare our re- 
sults with those of Kapusta, Lichard, and Seibert [[l3[ , 
which did not include the effects of the ai meson. 



g 







M 
\ 


a-2 






\ 


a=1.64 






a=1.05 


1 


_ 


.\ 


Kapusta et al _ 




- 




- 




- 


V 


- 




. 


\\ 


_ 




- 


% 


- 






\ Y' 


\ 






\ 


\\ 








X NX 


n 1 


_ 




'^ X. 












- 


7r7T->p7 


\ ^=5^..^ 












- 








- 




^ ^ ""^^ ~ 




- 




^\^ ^ 


ni 


- 


1 


"^ 



0.5 



1.5 



(GeV) 



FIG. 2. The total cross section for nn -^ p'y for a 
1.64 and 1.05. The value of g is determined by nip 
The result of Ref. ||l3] is shown for comparison. 
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For TTTT -^ p7, the a = 2 and a = 1.64 results are 
similar and near the results of Kapusta, et ai, within the 
pertinent values of ^/s. Likewise, the irp —^ ttj result 
for a = 2 (which has a visible bump close to ^/2nip) 
drops quickly to the Kapusta, et al, cross section due to 
the broadening in the momentum dependent ai width. 
The a = 1.64 result however, is larger by almost a factor 
of two, showing there is no simple connection between 
having a small Ta-^—t-yn width and the cross section in 
the np channel. Therefore we will use a = 1.64 for the 
free space calculations below and not concern ourselves 
further with the a = 2 case. As a consistency check, 
our results reduce essentially to those of Kapusta, et at, 
when we turn off the ai effects, as they should. 

The a = 1.05 results lead to substantially smaller cross 
sections, indicating that less photons will be produced 
from these processes when chiral symmetry restoration 
is approached. This is especially seen when comparing 
to the free space result (a — 1.64). Including the scal- 
ing of the parameter g as well will lead to a decrease of 
the vector meson masses as discussed in the next sec- 
tion. The suppression of cross sections seen here will be 
counterbalanced by an enhancement from the increase in 
phase space and the rates will not appreciably change. 
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FIG. 3. The total cross section for np — > 7r7 for a — 2, 
1.64 and 1.05. The value of g is determined by m^ = ag^f^. 
The result of Ref. ||l3] is shown for comparison. 

Emission rates for photons from a gas of hadrons at 
thermal equilibrium can be readily computed from the 
matrix element for a given process. For the process 1 -I- 
2^3-1-7, the differential rate is given by 



E. 



dR 
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(fpi (Pp2 



d^P3 



\M\' 



'^(Pp^ 2(27r)3 J {2Tif2Ei {2tt)^2E2 {2TTf2E-i 
i27r)^S\p,+p2-p3^P-,)f{Ei)fiE2) (1 + fiEs)) 

where A4 is the scattering amplitude, f{x) — 
(exp(a;/T) — 1)~^ is the Bose-Einstein distribution, and Af 
is the degeneracy factor which equals the number of dis- 
tinct incoming states considered in the matrix element. 
Since the matrix element depends only on the Mandel- 
stam variables s — {pi +P2)^ and t = {pi —p^Y', we follow 
]l3[ and insert integrals over each of these variables with 
delta functions ensuring the above definitions. The dif- 
ferential rate can then be expressed in terms of an exact 
quadruple integral as shown in Appendix B. 

For the thermal emission rate from the decay p -^ mrj, 
we used the approximate expression given in [ [L3[ . As it is 
seen from Figs. Q, ^ and ^, our thermal rates for a — 1.64 
are considerably different from those of Song ||l^ who 
included the oi meson within the massive Yang-Mills 
approach but did not include a momentum dependent 
width. For the tttt -^ p7 process, the rates are similar 
to Song for small photon energy but soon fall below. Al- 
though the massive Yang-Mills Lagrangian is a special 
case of HLS, we have very different parameter choices. 
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FIG. 4. Thermal rates for the process ivp -^ n"/ as com- 
pared to those of Song and Kapusta et al. 



width. At thermal equilibrium, the effects from this de- 
cay are outshined by those due to the two-body processes 
by almost two orders of magnitude. 



> 



•a 




10-10 
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The variable ai width is an important effect, leading 
to a very broad resonance in the s-channel irp —^ n^ cross 
section. There the contribution of the ai is practically 
zero for ^/s > 1.4 GeV. However, this does doesn't ac- 
count for the factor of approximately 3 between our result 
for the thermal emission rates and those of Song. For the 
p decay, our result is almost identical to Song's, as ex- 
pected from our agreement on the corresponding decay 



IV. SCENARIOS FOR THE CHIRAL PHASE 
TRANSITION 

The most attractive feature of HLS is that it describes 
fairly well the dynamics of pions and vector mesons with 
only three parameters: the pion decay constant /^, the 
universal coupling constant g, and the parameter a. From 
the construction of the HLS Lagrangian, the latter two 
define the masses of the vector mesons 



2 2 j-2 



ml/ml^ = (a- l)/c 



(3) 



In practice, this relation may be reversed in order to use 
the vector meson masses to set the values of a and g. 

We set out to use this description for the physics of 
mesons at high temperature and baryon density. Our 
basic assumption is that, at least for moderate condi- 
tions, the HLS description still holds, only with changed 
parameter values. In other words, the effects of tempera- 
ture and medium may be described in terms of changing 
U, rup, and ma^ only. 

Undoubtedly, if our description is valid, renormaliza- 
tion effects as well as temperature and medium effects 
within the HLS description pOl play a major part in the 
evolution of our three parameters. On the other hand, 
there may be more fundamental effects, which can not 
be predicted from within the effective theory. After all, 
spontaneous chiral symmetry breaking is a property of 
full QCD, and the mechanism controlling it may or may 



not be captured in an approximate description like ours. 
Therefore, we do not address within the HLS approach 
the issue of how the three quantities evolve towards the 
phase transition. Instead, we will consider a number of 
scenarios based on very general arguments. 

The p meson mass has been suggested as a possible or- 
der parameter for the chiral phase transition |g^ . If the 
phase transition is of second order, then nip will smoothly 
decrease towards zero. According to the scaling argu- 
ment by Brown and Rho ||8||, all other quantities such as 
/^ and TOai will be driven by the ratio m* /nip {m* is 
the in-medium p mass). Our aim is to find out the ef- 
fect of such a picture on photon emission from heavy ion 
collisions. 

First scenario: The conventional wisdom would be to 
allow the oi mass to drop towards zero along with the p 
mass 



nrip niai 

However, this simple scaling will freeze the HLS a param- 
eter at its free-space value as can be seen from eq. (||). 
This implies a will never attain its renormalization group 
fixed point of the HLS theory [Q , a = 1 , which would 
lead to the Georgi vector limitR. The dependence of g on 
the in-medium p mass is determined by the scaling di- 
mension of /^ as given by the Gell-Mann-Oakes-Renner 
relation: 



1/2q 






The suggested values for a are 1/3 |], 1/2 J2J, and 1 
p5| . Since g scales like (ipip)"'^^/^ from the above and 
eq. pi choosing a = 1/3 would lead to an infinitely strong 
coupling in the chiral limit contrary to the asymptotic 
freedom of QCD. The choice of a = 1/2 is therefore a 
limiting case, leading to 

/tt riip 

and thus g is frozen to its free-space value just like a is. 
This prescription is equivalent to the one used in [|l7| . 

Second scenario: As opposed to the complete disregard 
for the renormalization group behavior of a and g, we 
now take the fixed points to be realized. If a approaches 
1 as the p mass goes to zero, the ratio ma^lrtip must 
diverge, as it can be seen from (0). This does not imply 
the fli mass can not vanish, only that it must vanish 
slower than rup. As a limiting case, however, we take 



mai to be constant]^ and this determines a. Note that 
now the flow of a and g to their fixed point values of 
1 and respectively are intimately connected such that 
they conspire to give a constant rua^ as seen in eq. (^ . Of 
the possible values cited, only a = 1 for the Brown-Rho 
scaling drives g to 0. In conclusion, we take 



This is one way of "implementing" the Georgi vector 
limit, without having to address the issue of the fate of 
the a\. If the a\ becomes massless along with the p, then 
its longitudinal polarization will be a massless scalar, de- 
generate with the pion. 

Third scenario: A different physical picture that may 
be considered [g6| as an alternative to the vanishing of 
nip in the chiral limit is based on lattice results indi- 
cating that mp and nia^ decrease with temperature but 
become equal at some finite value, not running to zero 
p2[ . If a similar situation occurs with density, we may 
parameterize it by taking the vector meson masses rup 
and ma^ as the sum of an invariant piece nrideg, which is 
the same in both, and a piece that scales with the chiral 
order parameter to some positive power a: 



M - nideg V {ipip) 



for M = nip and m-ai • It is also possible to consider 
different values of a for the two vector mesons, but for 
simplicity we take them equal. Again, /^r will scale with 
0^/^" and so in the limit <f) —^ 0, g is proportional to 
0(a-i)/2a^ which imposcs the choice a = 1 among the 
ones mentioned above. Therefore /tt scales with 0^/^ and 
the matrix elements sharply increase as (p approaches 0. 
Although there are only a few occurrences of very low 
ni* they are significantly enhanced by this effect. In the 
following we take nideg = 0.2nip. 

Each of the three scenarios gives distinctly different 
results as illustrated in Figs. ^ and ||. Fig. |^ shows the 
dependence of the photon production cross sections (for 
fixed y/s — 1200 MeV ) and p radiative decay width on 
the p mass. Since both the strict Brown-Rho scaling of 
Scenario I and the lattice based scaling of Scenario III 
allow both vector meson masses to drop together, we see 
an increase in the cross sections as the p mass decreases. 
From this point of view, the main difference in the three 
scenarios is the way in which the oi mass decreases: ei- 
ther at the same rate as the p, not at all, or faster than 
the p. This leads to a different rate of enhancement for 



■^A one-loop calculation of the renormalization group equa- 
tions including the ai meson has not yet been done, but even 
a different fixed point a = oq will not change the final situa- 
tions described in this paper. 



■^In our simulations, we expect only relatively small effects 
from partial chiral symmetry restoration. Ifthemai decreases 
significantly slower than mp, we may take it constant in the 
first approximation. 



lower p masses, and no considerable change for the ex- 
treme case of the Georgi- vector limit in Scenario II. This 
can also be seen in Fig. || where the p mass is fixed at 
500 MeV. The ai mass reduction in Scenario I and III is 
accompanied by an increase in the cross sections. 
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energies. In a more careful analysis, one should interpo- 
late between the effective theory and asymptotically free 
QCD. For our purposes, we calculate our cross sections 
up to where they start to increase again {^/s = 2 GeV 
for the free space case) and then fix the cross section to 
that final value for larger values of y/s. In Fig.|8|, we in- 
dicate this by the thicker lines below the cutoff y^s. As 
long as the medium effects we are considering are rela- 
tively small, the cutoff value stays reasonably high, and 
our results should not be strongly influenced. However, 
one should have this in mind when looking at the high 
Pt region. 
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FIG. 8. Photon production cross sections in free space 
and for a p mass decreased to 500 MeV in the three different 
scenarios. 



FIG. 7. Photon production cross sections at ^/s = 1200 
MeV and p decay rate as a function of p mass, for the three 
different scenarios. 

Notice that in Fig. g the cross sections for reduced p 
mass tend to increase with y/s starting at approximately 
1300 MeV. The same behavior is present in the free-space 
result, only at higher values of s. This is certainly un- 
physical and comes from the terms with high powers of 
momenta in the extended HLS Lagrangian introduced by 
the redefinition of the A^ field. We are, after all, work- 
ing with an effective theory which breaks down at high 



Naturally, the thermal equilibrium rates shown at the 
end of Section III will increase dramatically for all three 
scenarios as the vector masses decrease. This is mostly 
due to the Boltzmann factor. For example, if the p-meson 
mass is reduced to 500 MeV at T = 150 MeV, the ther- 
mal emission rate pir —f 777 increases by about a factor 
of 6. However, as we will discuss in the next section, 
this effect is not the only factor to consider when ap- 
plying these results to heavy ion collisions. In addition, 
the total pion multiplicity is constrained by experiment. 
Whereas the thermal rate calculation corresponds to a 
grand-canonical ensemble of mesons which can produce 



excess pions at will, the real situation in heavy-ion col- 
lisions is closer to an ensemble with fixed pion number 
constrained by the hadronic observables. 



V. PHOTON SPECTRA IN HEAVY-ION 
COLLISIONS 

In studying medium effects in heavy-ion collisions, 
relativistic transport calculations [pTpq] based on the 
Walecka-type model has been quite useful, providing a 
thermodynamically consistent description of the medium 
effects through the scalar and vector fields. In heavy- ion 
collisions at CERN-SPS energies, many hadrons are pro- 
duced in the initial nucleon-nuclcon interactions. This is 
usually modeled by the fragmentation of strings, which 
are the chromoelectric flux-tubes excited from the inter- 
acting quarks. One successful model for taking into ac- 
count this non-equilibrium dynamics is the RQMD model 
p9| . To extend the relativistic transport model to heavy- 
ion collisions at these energies, we have used as initial 
conditions the hadron abundance and distributions ob- 
tained from the string fragmentation in RQMD. 

Further interactions and decays of these 'primary' 
hadrons are then taken into account through a con- 
ventional relativistic transport model. We include non- 
strange baryons with masses below 1.72 GeV, as well as 
A, A(1405), S, and E(1385). For mesons we include tt, 77, 
/9, uj, 77', ai, and cj), as well as K and K*{892). Baryons are 
propagated in their mean fields, which are assumed to be 
the same for all non-strange baryons. The mean fields for 
hyperons are assumed to be 2/3 of that for non-strange 
baryons, based on the simple quark counting rule. 

In addition to propagation in mean fields, hadrons 
also interact under stochastic two-body collisions. For 
baryon-baryon interactions, we include both elastic and 
inelastic scattering for the nucleon, A(1232), A^(1440), 
and A^(1535). Their cross sections are either taken from 
Refs. [ P0|,^ or obtained using the detailed balance pro- 
cedure |32P The meson-baryon interactions are modeled 
by baryon resonance formation and decay. For exam- 
ple, the interaction of a pion with a nucleon proceeds 
through the formation and decay of any of the N or A 
resonances from the A(1232) up to the iV(1720). The 
formation cross sections are taken to be of the relativis- 
tic Breit-Wigner form. The meson-meson interactions 
are either formulated through resonance formation and 
decay when the intermediate meson is explicitly included 
in our model, such as the ai meson, or treated as a di- 
rect elastic scattering with a cross section estimated from 
various theoretical models [g[. 

Photon creation is taken into account during the evo- 
lution of the transport code through decays or two body 
processes. The experimental data for lu and ry' decays 
are used as these particles are long-lived and mostly de- 
cay after escaping into free-space. Photons can also be 
produced from the decay of baryon resonances. These 



contributions are usually neglected in hydrodynamical 
calculations |3^,Q, but are included here through ex- 
perimentally measured radiative decay widths pl| ] . 

Otherwise, we include the three main contributors to 
photon production: p —^ tttt"/, tttt —^ frf, and irp —^ 717 
as calculated in the previous sections and modify the two 
fundamental parameters a and g of the HLS Lagrangian 
according to each of the three scenarios mentioned above. 

In Fig. H we compare our thermal photon spectra with 
the upper bound of the WA80 collaboration. Overall the 
results are below the upper bound for photons with trans- 
verse momenta well below 1 GeV. For higher transverse 
momenta, our results touch the upper bound of the ex- 
perimental data. These observations are essentially the 
same as those of Ref . |Q , although the contribution from 
two-body collisions at high transverse momenta now be- 
comes comparable to that of meson decay, due to a more 
realistic treatment of oi effects. This shows that whether 
the vector meson masses are reduced or not, the rates do 
not change much since the opening up of phase space is 
balanced against a decrease in the coupling constants. 
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FIG. 9. Total thermal photon spectra in central S-|-Au 
collisions at 200 AGeV for the different scenarios discussed. 

The fact that there is no dramatic increase in pho- 
ton emission rates would seem to contradict the naive 
expectations of thermal rate calculations. Regardless of 
whether the dropping mass scenario is invoked or not, 
any photon spectrum calculation must at the same time 
fit the observable pion multiplicity. Through the use 
of a complete transport calculation as described above, 
this quantity can be consistently obtained. As shown in 
Ref. |g] , this could be achieved in a thermal model by ei- 
ther decreasing the vector meson mass or increasing the 
pion chemical potential by hand. Thus, for a meaningful 



comparison with the dropping mass scenarios, one would 
need to use a large chemical potential in the bare mass 
scenario, which would push up the thermal rates. In this 
paper however, we do not use thermal equilibrium ini- 
tial abundances. Instead, we evolve the system from the 
same initial RQMD output for all cases considered. This 
procedure guarantees the correct final pion yield without 
additional assumptions such as a pion chemical potential 

With a dropping p mass, the p number increases in the 
initial stage mainly through the process p <-*■ tttt. This 
however is also the dominant decay mode for the p and 
so the pious are eventually regained at freeze-out. The 
three pmrj photon-producing processes discussed in de- 
tail above are not enhanced dramatically if the number of 
p mesons increase because this is achieved at the expense 
of reducing the number of pions. As a result, pn — > t:^ 
in neither favored nor disfavored and the other two pro- 
cesses TTTT — > pj and p — *■ 7r7r7 are balanced against each 
other. Although all three scenarios are close to the free 
space result, the largest change comes to Scenario II and 
the Georgi vector limit, which is reduced by almost a 
factor of two in the final rates as seen in Fig. O. 



1. The lowest rates come from Scenario II as one would 
expect from the cross-sections. 
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FIG. 10. The same as Fig. 0, for p -^ tv-kj. 

This variation in the three scenarios can be seen bet- 
ter by breaking the rates up into the important contri- 
butions. This is shown in Figs. |l^, |Tl|, and |2| for the 
p —f 77777, TTTT —^ p7, aud TTp — > 777 ratcs respectively. 
The main contribution can be seen to come from p de- 
cay, which is larger than what one would expect from 
the thermal rates of Section III due to off-shell propaga- 
tion enhancement. The largest value to the decay rates 
comes from Scenario III in which both masses drop, but 
a increases as the ratio of the p and ai mass approaches 
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FIG. 12. The same as Fig. g, for yrp — > 777. 

It is also interesting to note how similar Scenario I is to 
the free rates in all three figures. This is because Brown- 
Rho scaling for the ai meson does not allow either of the 
two parameters of the HLS Lagrangian to change. The 
widest range in values among the scenarios comes in the 
7777 — !■ /97 rates which give a spread of almost an order of 



magnitude between Scenario III (largest) and Scenario II 
(smallest). One can see from Fig. ||that this process has 
a larger spread between scenarios in terms cross sections, 
too. 

It is interesting that in total, all three scenarios as well 
as the free rates touch the upper bounds set by WA80 
for large pt- If other processes such as oitt -^ p7 were 
included, they would also feed into this high momentum 
region and possibly push the rates past the upper bounds. 
One must be cautious however in drawing conclusions in 
terms of absolute values of the photon rates at high pt 
even though we cut off the unphysical high-momentum 
dependence of our cross sections as discussed at the end of 
Section IV. Needless to say, results from WA98 and future 
photon measurements are needed to properly interpret 
the discrepancy with the present data. 

With excellent data in the photon spectra, one could 
possibly distinguish between the three broad scenarios set 
forth in this paper. However, given the extremely poor 
signal to noise ratio for the photons, the experimental 
uncertainty is too large to say whether one or another of 
the curves in Fig. O is right. One might rank the scenarios 
in order of decreasing total photon rate. Then, the free- 
space result and Scenario III seem more likely to exceed 
the upper bounds for the region pt < 1.75 GeV, than 
Scenarios I and II. Both the free-space result and Scenario 
III are strong in this region for p decay. The free-space 
result outshines all other scenarios for pir -^ 777, while 
Scenario III dominates for tttt — > p^. Except for p decay 
above 2 GeV, Scenario II gives the lowest predictions. 
Also, it is the only one to have the same slope as the 
WA80 hmits. 



VI. SUMMARY AND CONCLUSIONS 

The purpose of this study is to investigate the effect 
of dropping in-medium vector meson masses on photon 
spectra from heavy ion collisions. This has been done 
earlier for dilepton spectra H and for photon spectra 
|l7| . In the case of photon spectra, the effect of the ai 
has to be included p4[, since it is potentially important. 

We recalculated three processes important for photon 
production in heavy ion collisions using the extended hid- 
den local symmetry Lagrangian |18[ with the proper oi 
mass. We find that the contribution of the ai as inter- 
mediate state is relatively small, in contrast to results 
found in earlier models plJl^, but closer to results in 
Ref. |l6J. It would be interesting to see whether this is 
indeed a genuine physical effect by looking at alternative 
formulations of HLS. In particular, the ir-ai mixing is 
taken care of now by a shift of the oi field whereas using 
one of the unphysical fields that appear in non-unitary 
gauges could have the same effect without modifying the 
aiTTj vertex. Both the ai —^ nj decay width and the s- 
to d-wave ratio in ai decay should be used to identify the 
most advantageous approach. 



However, the small effect of the ai meson lies in some- 
thing even more basic: the linear relation in the HLS 
between the coupling g and the vector meson masses. 
Therefore, as temperature and density increase, the 
masses drop in agreement with Brown-Rho scaling |q| 
and the couplings are driven to zero in accordance with 
asymptotic scaling. This suppresses the rates when den- 
sity corrections are taken into account. 

Nevertheless, using only the vector and axial- vector 
meson masses as input parameters our model reasonably 
predicts the vector and axial- vector meson decay widths. 
Our predictions for thermal photon emission rates from a 
hadronic gas through the processes pir — *■ ttj, tttt — > p7, 
and p -^ TTTTj are within the range of similar results in 
the literature. The rates for pn —f nj are slightly higher 
than those of Kapusta, et al. [0| , but a factor 2-3 lower 
than those of Song |g5| . 

We considered three different scenarios for the in- 
medium evolution of nip and ma^ in an attempt to un- 
derstand the differences in signals they would give: strict 
Brown-Rho scaling where the ratio nip/mai stays fixed, 
an extreme case of the Georgi- vector limit where niai is 
fixed, and a schematic model inspired by finite tempera- 
ture lattice results. 

We simulated single photon spectra in central S+Au 
collisions at SPS energies using the relativistic transport 
model that has been used to study dilepton spectra in 
the same reactions. We included photons from the back- 
ground sources of tt'^ and 77 decays, as well as thermal 
sources such as meson decays, decays of baryon reso- 
nances, and two-body processes. We found that more 
than 95% of single photons come from the decays of 7r° 
and 77. The thermal photons account for only less than 
5% of all single photons, in agreement with the exper- 
imental observation made by the WA80 collaboration. 
We compared our thermal photon spectra with the ex- 
perimental upper bound extracted by the WA80 collab- 
oration. Overall, the results were all comparable to each 
other for the total rate, none of them exceeding signif- 
icantly the experimental bounds. For pt < 2 GeV, the 
largest photon yields were found for the lattice-inspired 
scenario and the simulation without dropping masses. 
The Georgi- vector limit case gives the lowest yield in this 
range. It slightly exceeds the upper bounds of WA80 for 
large pt, but a fully consistent high energy behavior is a 
matter of further study. 

In conclusion, we find that in the extended hidden local 
symmetry approach the role of the oi for the processes we 
considered is less important then it is generally thought. 

We showed that one can implement the dropping in- 
medium p mass and Brown-Rho scaling without violating 
the existing experimental limits for photon production. 
Furthermore, there are indications that different scenar- 
ios lead to experimentally distinguishable predictions. 
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APPENDIX A: MATRIX ELEMENTS 

For p'^{p) -^ TT^{pi)TT''{p2)j{k), the matrix element A4 shown diagrammatically in Fig. n^is given by the addition of 
the following four expressions (r — ml/fri^^). 

Ml = e"''"'e^'"'--—- {m^ei-e2 + 2r {ei-e2P2-p - e.i-P2(^2-p) + 2r^ {ei-e2Pi-p2 ~ ei-Pie2-P2)} + (pi,&) ^ {P2,c) 

a3e bee e j , , , , . ^API + P2)^ + mf, 1 , / i,N / N 

M2 = €" e^^—p: \ [ei-pie2-p + ei-ke2-pi - ei-e2Pi-k) — • f- - ^r^i -Pi e2-P2 > + {pi,b) ^ {p2,c) 

gj4 I [Pi +P2Y ~m^p I 



,aee,3&e Cj -pa £2 'Pi 

(pi + fc)2 - ml 



M3 = 2age (1 + C.) ^"''^'" ^^ ,,.2- 2 + (^i' ^) ^ (^2, c) 



_ ipi-ke2-W-e2-pik-W) \ ^ 
ipi + kf - ml^ J 



with ^r = 2r(l — 2r) and 

W^^ - ei -^2 [(4r - IK - rP2 ] - [(4^ - 1)P2 -p + rpa • (pi + fc)] e^'. 
As required by gauge-invariancc, the total matrix clement M vanishes under the replacement £2 —^ k. 

APPENDIX B: DIFFERENTIAL THERMAL RATE 

For the process 1 + 2^3 + 7 we define s = (pi +P2)^, t = {pi — p^Y . With E3 = Ei + E2 — E^, introduce 

qi^ml-t; q2 ^ s + 1 - ml - mj; q^ ^ s - mj ^ qi + q2; p\ = -ml + -—qu - ql, /c = 1,2,3. 

The differential rate is then 

<!« _ ^' f,. f,., .,,-.., 2 f,^ [,^ S {El) S {£■,){! + 1(E,)) 



^^^ ^ i^^nus, h h |^<"' '*!' Z"^' hrf. 



A i2 ,2 i /2 /2 /2\^ ' 

4p'iP'2 - (P3-P1 -P2) 

The integration limits on s and i, in addition to s > (ttt-i + 777,2)^, are such that gi > and 52 > 0. The limits on i?i, 
E2 are set by 

_ q\ Em^ „ (72 Em'^ ^ ^ „ 

Ei> ^ + ^; ^2 > tI- + ^; Ei+E2> E^. 
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